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£N| , Abstract. A relatively recent study by Mars and Senovilla provided us with 

i— h ' a uniqueness result for the exterior vacuum gravitational field generated by an 

i isolated distribution of matter in axial rotation in equilibrium in General Relativity. 

The generalisation to exterior electrovacuum gravitational fields, to include charged 
—i ■ rotating objects, is presented here. 
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^ ' 1. Introduction 

o : 

in 

A proper understanding of rotating objects in equilibrium within the context of General 
Relativity is fundamental for many astrophysical situations. Very unfortunately, finding 
global models for rotating objects in General Relativity has proven to be extremely 
difficult, even for axially symmetric configurations in equilibrium. So far, there 
are no known complete explicit models for a self-gravitating finite body (with non 
^ ■ empty interior) together with its exterior except for spherical, and hence non-rotating, 
configurations. 

This work focuses on the theoretical point of view of the construction of global 
models for finite objects by means of the matching of spacetimes: the whole configuration 
is composed of two regions, one region of spacetime (W 1 ,^ 1 ) describing the interior (I) 
of the body and another (W E ,g E ) describing the exterior field (E), matched across 
a hypersurface E. This matching hypersurface is then a common boundary of the two 
regions and corresponds to the limiting surface of the body at all times. The two regions 
can then be treated independently, taking into account that the two problems will have 
to satisfy compatible boundary conditions on £ imposed by the matching conditions. To 
model the equilibrium state of the rotating configuration, the whole matched spacetime, 
and hence each of the interior and exterior regions, is taken to be strictly stationary. In 
addition, it has usually been assumed that the model is axially symmetric. To account 
for the isolation of the body, the exterior region is required to be asymptotically flat. 

In previous works on models of isolated rotating bodies in equilibrium, only vacuum 
exteriors have been taken into account. Nevertheless, compact objects in astrophysics 
have been also considered many times as sources of electromagnetic fields (see and [3] 
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and references therein). In fact, the existence of a net charge, no matter how negligible, 
has drastic consequences for the global spacetime structure, in principle. The aim of 
the present work is to revisit the exterior problem, and generalise the uniqueness results 
found for the vacuum case in [2] to include stationary and axisymmetric electromagnetic 
fields without sources (electrovacuum fields) . The main result deals with the uniqueness 
of the exterior electrovacuum field generated by an isolated distribution of charged 
matter in axial rotation and in equilibrium. It is worth noticing that spherically 
symmetric global models for stellar collapse to include charged radiating stars and voids 
has been already considered by several authors during the past three decades (see |3] 
and jlj and references therein). 

The exterior stationary and axisymmetric Einstein-Maxwell electrovacuum problem 
consists of an elliptic PDE system, known as the Ernst equations, for a couple of complex 
functions 0IE]; the electromagnetic potential A and the so-called Ernst potential £, 
which carry all the information about the exterior geometry and electromagnetic field. 
It is known that the matching of a vacuum exterior with a given interior, prescribing 
the identification of the interior and the exterior across S if required, determines: (i) 
the existence of the matching hypersurface, and if so, its form as seen from the interior, 
generally, (ii) the matching hypersurface as seen from the exterior, and (iii) the values 
of £ up to an arbitrary additive imaginary constant, say iu, and the normal derivative 
of £ there [21 Cj. This data constitutes a set of boundary conditions of Cauchy type for 
£ on S that depend on the parameter uj. 

In order to accommodate an electrovacuum exterior, it is necessary to, first, 
generalise the spacetime matching conditions for a vacuum exterior, also taking into 
acount the matching conditions for the electromagnetic field |Hj (see also jlj). For an 
object without surface charges^ the electromagnetic field has to be continuous across 
E, and hence the interior quantities determine the electromagnetic potential A up to 
an additive complex constant A. Here it will be shown how, when considering an 
electrovacuum exterior, the sets (i) and (ii) remain formally unchanged whereas the 
rest of matching conditions constitute a generalised new set (iii) of Cauchy boundary 
data for £ on £ in terms of the interior quantities and involving uj and A. Summing up, 
given the interior, the matching conditions determine Cauchy boundary data for £ and 
A on S up to three degrees of freedom, given by uj and A. 

Due to the elliptic character of the exterior electrovacuum problem, the Cauchy 
boundary data, together with the decaying conditions at infinity, constitute an 
overdetermined set of boundary conditions for the exterior electrovacuum problem. 
Nevertheless, as shown later, the three parameters encoded in u, A are in principle 
physically relevant, and therefore, the exterior problem is overdetermined but, a priori, 
not unique. The proof of the uniqueness of the exterior solution has to follow two steps, 
then. In a first step it is shown that given the values of £ and A on £ (Dirichlet boundary 
data), the exterior field is unique. On a second step, it is shown that the completion to 

| This assumption is not really necessary, since one could simply consider a surface charge and then 
include that as information that has to be given with the rest of the interior configuration. 
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a whole set of Cauchy boundary data of a Dirichlet boundary data set that depends on 
lu, A determines the values of uj and A, and hence the Dirichlet data, provided that the 
solution exists. 

The paper is structured as follows. Section El is devoted to a brief review concerning 
the electrovacuum exterior problem. The dependence on the free parameters u and A 
of the Cauchy boundary conditions on £ is determined in Section El For the sake of 
brevity in the main text, the full determination of the matching conditions is left to 
|Appendix A| Section H] deals then with the proof of the uniqueness for the Dirichlet 
problem. Finally, Section El is devoted to showing that the full set of Cauchy boundary 
data determine the values of uj and A if the solution exists. The statement of the final 
result is left to a conclusions section. 

The units all throughout the paper are chosen so that G = c = 1. 

2. The exterior problem 

For completeness and to fix some notation, this section is devoted to a brief review of 
the stationary and axisymmetric electrovacuum problem. Let me refer to e.g. [SI Ej for 
the details. 

Given the electromagnetic field described by the 2-form F, it will be convenient to 
use the self dual 2-form F defined as 

J= = F + i*F, 

where * stands for the Hodge dual. The self-dual 2-form J- is then completely 
determined by an arbitrary non-null vector k and the projection f ^ = J-(-,k,), which 
clearly obeys f g (n) = 0, by 



where we have defined = k- 0), and A denotes the exterior product (see e.g. [Hj). 
The real and imaginary parts of the form fg are usually denoted by Eg and Bg, and 
for unit timelike k, they correspond to the electric and magnetic parts of F with respect 
to k, respectively, i.e. [3] 



where j is the electromagnetic current source, and the energy-momentum reads 



where the bar denotes the complex conjugate. 

Let us consider a strictly stationary§ and axisymmetric [3] asymptotically flat, 
causal and simply connected C 3 spacetime with connected boundary (W E , g E , S E ), 

§ In global terms, we demand that the manifold admits an everywhere timelike Killing vector with 
complete orbits. 



^=Tr[/*AK + i*(/ g AK)] 



fg = Ex + iBit. 
The Maxwell equations in terms of 7- take the form 




(1) 
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containing an electrovacuum field J- invariant under the same group of isometries, 
i.e. the stationary and axisymmetric pair (g E ,^F) is a solution of the Einstein-Maxwell 
equations without sources (j = 0). The intrinsically defined axial Killing vector field |Hj 
(see also 0) will be denoted by ff, and the axis of symmetry by W^, which is non-empty 
by assumption. £ will denote an everywhere timelike Killing vector field which together 
with ff generate the necessarily Abelian ^U] (see also [3]) G% group of isometries that 
act on timelike surfaces T 2 . It is also well known that the G<i on T 2 group of isometries 
in a stationary and axisymmetric electrovacuum spacetime must act orthogonally 
transitively (see e.g. j3j). These properties, together with the vanishing of the trace 
of the energy-momentum tensor for the electromagnetic field, imply the existence of a 
coordinate system {t,<f), p, z} with ranges within p > 0, t, z e R, < 4> < 2n, the so- 
called Weyl canonical coordinates, in which the line-element for W E \ W% reads locally 

B 

ds 2 = -e 2U (dt + Ad<P) 2 + e~ 2U [e 2k (dp 2 + dz 2 ) + p 2 d<p 2 ] , (2) 

where U, A, k are functions of p and z, the axial Killing vector is given by ff = d<p, 
the axis of symmetry W% is located at the limit p —>■ 0, and £ = d t . Moreover, the 
coordinate t can be chosen to measure proper time of an observer at infinity. In that 
case the Killing vector £ is intrinsically defined by being unit at infinity 5\. Within this 
setting, the coordinate freedom in (J2J) consists only on trivial constant shifts of t, <fi and 
z. 

In order to have our G 2 group defining complete orbits, the Killing vectors {£, if} 
are required to be tangent to S E . In other words, we demand that E B preserves the 
stationarity and axial symmetries ^T]. As a consequence, E E is everywhere timelike by 
assumption. Furthermore, since W E is assumed simply connected, connectedness of S E 
implies that the slices of constant t of S E are homeomorphic to two-spheres. 

The existence of the Weyl canonical coordinates globally on W E \ W% had been 
ensured in relation with the black hole uniqueness theorems ^2] (see also pU), where 
S E lies on the horizon, so that p|sE = 0. In the present case, though, the boundary 
S E lies on the surface of the object, and therefore those results cannot be used in 
principle. Nevertheless, the good behaviour of the function p can be also ensured when 
the boundary S E at constant t is an axially symmetric surface homeomorphic to a two- 
sphere, as it is the case here. Global Weyl coordinates had been used in previous works 
on global models of finite objects, but without proof 0IZI- Due to its length, and not to 
overwhelm this paper, it has been preferable to leave the proof on the existence of the 
global Weyl coordinates in the present scenario to a separate work [Hj. Taking global 
Weyl coordinates (J2J) for (W E , g E , S E ), the axially symmetric hypersurface S E is thus 
given parametrically by 13 Clj 

S E : {t = T,<p = <p,p = p(ji),z = z(fj)}. (3) 

Denoting by k any of the two Killing vectors £ or ff, the invariance of the 
electromagnetic field F with respect to k, equivalent to L^T- = 0, implies the existence 
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of a complex scalar potential Ag such that (see e.g. 010]) 
dA« = /« 

for either k. Moreover, for stationary and axisymmetric electrovacuum fields, f u must 
be orthogonal to the orbits, and thus = A^(p, z) E]. 

The Einstein-Maxwell equations free of sources in the simply connected manifold 
W E imply, in turn, the existence of a real scalar £1%, the so-called twist potential, such 
that [inj (see e.g. jS]) 

dA* = Wjt-i (AgdA* - A^dA*) , (4) 

where tug e*(kA d/c) is the twist 1-form of Recalling the definition = (k • k), 
the so-called Ernst potential £ % (with respect to k) reads [TOJ El U2] 

£k ATg A^A^ -|- 

The Einstein-Maxwell equations now reduce to the elliptic system of equations for the 
potentials kg and £%, known as the Ernst equations, plus a quadrature for the remaining 
function in the metric (k(p, z) in (j2J) in terms of the potentials (see |Appendix A| for the 
explicit expressions). Dropping the k subindexes, the Ernst equations read JSJE] 

N5 ij di(pd 3 £) + p 5 ij di£ (d 3 £ + 2kd 3 k) = 0, (5) 

N5 ij di(pd 3 k) + p 5 i] dik (d 3 £ + 2A<9 i A) = 0, (6) 

where N = —{£ + £ + 2AA)/2, indices z, j correspond to {p, z}, and S i3 represents the 
2x2 identity. The solutions of the Ernst equations for k = £ and k = ff, that is A^*) 
and (£fj,kff), are called conjugated to one another and are bi- uniquely related (see |G]). 

Given either Ag) solution of the Ernst equations all the information of the 
exterior electrovacuum solution is recovered. In particular, the metric for (W E , g E ) in 
the form (J2J) is directly obtained by taking k — £ (= d t ): U is obtained from 

N^=-e 2U , (7) 

and A is determined, up to a constant, by the quadrature 

where * is used to denote the particular Hodge dual on the {p, z} two-plane, so that 
-kdz = —dp, *dp = dz, and is in turn obtained from taking the imaginary part of 
d£g and using (j3J). The function k is then fixed, up to an additive constant by the 
aforementioned quadratures involving £? and k^ j^j. 

The boundary conditions for the exterior problem consists of the conditions on the 
boundary associated to the surface of the body plus decaying conditions at infinity. 
Asymptotic flatness on (W E ,g E ) determines the behaviour of the Ernst potential at 
infinity ^7j (see also PHI Ej and references therein) and, in adition, requires asymptotic 
conditions on the electromagnetic field jTTj. The conditions at infinity read 

Ug = -Mr' 1 + 0(r- 2 ), ^- = -2zJr^ + 0(r" 3 ), (9) 
k^-Qr^ + Oir- 2 ), (10) 
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where r = yj p 2 + z 2 , for some constants M, J and Q. 

The next section is devoted to the boundary conditions on the hypersurface E that 
result from the matching conditions with a given stationary and axisymmetric interior. 

3. The matching conditions: boundary conditions on E 

Let us consider a given stationary and axisymmetric spacetime with boundary 
(W 1 ,*? 1 , E 1 ) containing an electromagnetic field ZF X , and describing the interior region 
matched across the stationary and axisymmetric (and thus timelike) hypersurface 
E(= E 1 = E E ) to the stationary and axisymmetric asymptotically flat electrovacuum 
spacetime (W E ,g E , E E ). Since we do not consider any superficial charge on the surface 
of the body, the usual junction conditions of electrodynamics jS] imply that the 
electromagnetic field has to be continuous across the timelike hypersurface S, and thus 
(see e.g. @]) 

^ I |s = ^ r |s. 

This clearly leads to a relation for the gradient of the electromagnetic potential at the 
exterior with respect to any vector field k, A#, and the interior quantities, given by 

Ag| s = e 3 a ^> /3 | E , ™( A *)ls = n a r a ^\ s , (11) 

where the dot denotes differentiation with respect to p, e 3 denotes the corresponding 
vector field, tangent to E, and n is normal (not necessarily unit) to E. It must be stressed 
that the absence of a surface charge is made here only for simplicity reasons. One can 
always consider a surface charge and simply add it to the information required from the 
interior region. In any case, the data would then be fixed by the electromagnetic 
matching conditions (see j3j). In short, given the interior electromagnetic field, the 
matching fixes 

dA^| s , (12) 

for either k for the exterior problem. As a consequence, Ag|s are fixed up to 
transformations of the form 

A^| s - A'a| E = A*| E + A*, (13) 

where is an arbitrary complex constant for each respective k. 

Regarding the matching of the spacetime (W E ,g E ) with a given (W 1 ,^ 1 ), it is 
proven in |Appendix A| (see also [21 E| for particular cases) that the matching conditions 
determine the existence and thereby the form of the matching hypersurface E = E 1 = E E 
as seen from both sides, in the general case, together with the values of the metric 
functions U and A (see (j2J)) and their normal derivatives on E. It is worth noticing that 
the additive constants in A and k (see above) are then determined. Also, one still has to 
keep in mind that it has been assumed that the identification of the exterior and interior 
across E has been prescribed, in order to fix two extra degrees of freedom introduced by 
the matching procedure, that correspond to the identification on E of the Killing vector 
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d t with a linear combination of two Killings vectors in the interior, see Appendix A and 



[3 EH d- The data on S for the exterior problem introduced by the matching with a 
given interior is thus given by 

Uk, dU\ s , (14) 
A| E , (L4| s , (15) 

which constitute Cauchy boundary data for the exterior problem. 

It is now straightforward to show how the data given by (|14 )) -([15 )) and (JT5)l translates 
onto data for 8. To begin with, the data (|14 )) -([15 )) implies that, given the interior, the 
matching conditions fix 

N x \ s , dN x \ s (16) 

for either k. Let us now consider the case k = £. Using the data (fTKJ) and (fTKj) . the 
relation (jHJ) implies that the interior geometry fixes Therefore, from and 

taking into account the freedom generated by A as introduced in ()1HJ) . the matching 
conditions fix dfi^s up to transformations of the form 

dtofc -> dfi£| E = d^| s - z (A r dA r - A r dA f ) | s , (17) 
and thus fix f^| s up to 

%|s -> = %|s + ^ - i(A { -A r - A r A € -) | E , (18) 
where cj|* is a real arbitrary constant. Although for the sake of clarity these relations have 
been explicitly deduced here for k — £, it can be shown that the matching conditions 
fix, in fact, 

With (19) 

for both k, and thus (|T7j) - (fTHj) hold for both k. Combining all the above, the Cauchy 
data on S for the Ernst potential is fixed up to transformations of the form (dropping 
the k subindexes) 

£|s = ^|s-AA + icu-2AA| E , and (20) 

d^| s ->d£% = d^| s -2AdA| s . (21) 

Due to the elliptic character of the Ernst equations system, the Cauchy data induced 
by the matching conditions overdetermine the problem. Nevertheless, as shown above, 
there are still three degrees of freedom on that data, generated by A and u. As remarked 
in [2], the parameter u is relevant, in principle, for the exterior problem, because the 
freedom in Q has been already used to set its asymptotic behaviour 0. For the same 
reason, expression (fTHj) makes clear that the parameter A is also relevant. All this means 
that although the exterior problem for a given interior is an overdetermined problem, 
the problem is still not unique. Therefore, in order to solve the uniqueness problem we 
have to address the following two questions: 

(a) uniqueness of the exterior solution {£ , A) given Dirichlet data on S, i.e. A| s }, 
and 
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(b) if a solution (£,A) for given data A| s , d£|s, dA| s }(A, u), fixed up to the 

transformations ()20 )) -()21 )) . exists, do A and u get determined? In other words, 
does the Cauchy data determine A and u, thereby fixing the Dirichlet data, if the 
solution exists? 

Clearly, an affirmative answer to both questions proves the uniqueness of the 
electro vacuum exterior field given a charged finite source. 

Regarding question (a), the procedure commonly used for proving uniqueness of 
Dirichlet problems is based on interpreting the Ernst equations as the Euler-Lagrange 
equations for a certain harmonic map. The proofs rely then, in particular, on the 
positivity of the Dirichlet functional associated to the harmonic map corresponding to 
the Ernst potential. This positivity is ensured for both k in the vacuum case, and the 
proof presented in j2] regarding question (a) is based on choosing k = £. However, in the 
presence of Maxwell fields (see below), one needs choosing k — ff. Adapting the proofs 
of electrovacuum black hole uniqueness theorems O E] to the present setting, Section 
E]is devoted to solve the uniqueness of the present Dirichlet problem using the Mazur 
identity approach. The choice of the Mazur identity approach over the Bunting identity 
approach used in [2j has been made for convenience, because the Mazur construction 
presented in Section |U is useful in dealing with question (b) in Section EJ 

The affirmative answer to question (b) acounts to showing that if a solution 
(£, A) for given data {£\s, A|s, d£\s, dA|g} exists, then another solution (£', A') with 
{£%,A%,d£'\x,dA%} given by dA'| s = dA| s , 0, (|2n jl -(|2T |l in terms of A and u 
exists only if A = = uj. This is done in Section Since this proof does not need any 
preference of vector k, I have preferred to present the calculations by choosing k = £, so 
that the corresponding quantities that appear have a usual and direct physical meaning. 

4. Uniqueness given Dirichlet data on X 

The proof follows those used in the uniqueness theorems of black holes (see ^21 E]), 
making use of the very rich intrinsic structure of the Ernst equations, the only difference 
being that the horizon is replaced by the matching hypersurface E. In what follows, the 
k indices will be omitted for simplicity in the expressions that hold for both k, whenever 
this does not lead to confussion. 

Equations ©-© can be interpreted as the Euler-Lagrange equations for the action 

m 

s = 4 / (i^\ de + 2AdA i 2 + ^i dA rt V, 

J w e \4N 2 * 1 1 / 

where \0\ 2 = g a ^6 a G/3 for any 1-form 0, and r\ is the volume element. Taking $ to be a 
hermitian matrix in SU(2, 1) defined by [0] (a, b : 1, 2, 3) 

$a& = Vab + 2 sign(A0^ 6 (22) 
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where v a = (2^\N\)~ 1 (£ — 1,2A, £ + 1) and r\ ah = diag(— 1, 1, 1), one can define a 
su(2, l)-valued 1-form by|| 

J = ■ d$. (23) 

The information carried by the solutions (£, A) has been now translated onto $. In 
terms of J ', the above action is rewritten as [T2J EH E| 



/w E 2' 

for which the variational equation reads ^21 E3 E] 

V Q X = 0. (24) 

Note that the way J7" is defined is not unique. In fact, although the definition given 
above corresponds to that in [12] after a trivial interchange J2 <-> J73 (see also P), later 
it will be more useful to refer to the construction presented in [TU] . 

Since the object $ is stationary and axisymmetric by definition, the 1-form J~ 
depends only on p and z and furthermore, it is "tangent" to the surfaces orthogonal 
to the orbits T 2 , i.e. J a = (0,0,^). Both $ and J" can therefore be taken as defined 
on the Riemannian space T that corresponds to the surfaces orthogonal to the orbits 
in (W E ,g E ). Recalling that S must be spatially homeomorphic to a two-sphere, the 
boundary of the Riemannian domain V of the solutions for $ is given by (see |14j): 

(i) A connected curve Tg, that corresponds to the projection of £ onto T, with ends 
at the axis, 

(ii) the two segments that correspond to either (connected) part of the axis, say T + 
and T_, where p = 0, 

so that dT = Ts U T_ U T + is connected. We will also use to denote infinity, at 
r 2 = p 2 + z 2 — > 00, this is, the ideal point of Y. Taking the metric 7 = \N\gf-dx l ® dx J 
on T one has V a J a = \N\p^ 1 ^V l (pJ'i), where denotes the covariant derivative 
with respect to 7, and thus equation reads 

(7) V 4 (P^) = (25) 

on (r, 7). The Dirichlet boundary data on Ts, i.e. $|x E , is obtained by construction 
from the data {£|x E , A| Te } which naturally corresponds to {£\s, ^|s}- 

The first key property for the proofs of the uniqueness theorems is the posit ivity 
of $ (and thus, of the action), which is ensured by taking k — ff (since sign(A^) = 1 
outside the axis) see ^21 E]- All the quantities in the remaining of this section will 
be associated to k = ff. Let us define now X = Nff (which in Weyl coordinates reads 
e~ 2U p 2 - e 2U A 2 > outside the axis), Y = and Aff = -Ae^X' 1 . The line-element 
(J2J) can be cast as 

ds 2 = -X' l p 2 dt 2 + X(d<f) + A^dt) 2 + X- l e 2k {dp 2 + dz 2 ). 
j Usual matrix product is denoted by {A ■ B) a b — J2 C A ac B c b and the trace reads Tr A = Y^, a A aa . 
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The problem with this choice is that X, which appears explicitly in the action, vanishes 
on the axis (located at p — 0), and therefore a careful analysis there is needed (see 
[T2*]). It is convenient to change the coordinates {p, z} to prolate spheroidal coordinates 
{x, y}, x > 1, \y\ < 1 by p 2 = v 2 {x 2 — 1)(1 — y 2 ),z = uxy, where v is an arbitrary 
positive constant, so that dp 2 + dz 2 = v 2 {x 2 — y 2 )ds 2 where 

1 1 

d 4 = ~n — 7 da; 2 + ^ dy 2 = 'j ij dx l dx :! . 

' x z — 1 i — y 

Taking the metric 7 = e~ 2k v~ 2 [x 2 — y 2 )~ 2 ^ on T, due to the conformal invariance of 
the equation for C (|23j). in (r,7) one has 

V\pJi) = 0, 

where "tilded" quantities will refer to the metric 7. In the {x, y} coordinates the 
boundary dT is given by T± : {y = ±1}, while is a connected curve for which 
x > 1 necessarily that joins y = 1 with y = —1, and one has : {x — > 00}. 

The core of the proof consists of considering two different sets of solutions $(1) 
and $( 2 ), with corresponding pairs , 3" (2) and field variables X^,X(2), etc, with 
common Dirichlet data on Ts, i.e. 

$(i)|r s = $(2)|t e - (26) 
Defining the differences / = — /( 2 ), for any functional / of the variables in $, and 

* = $(1) • $(2)~ 1 - I, 

where H is the identity, so that \& = -v4> $(i) = $(2), the Mazur identity follows (see 
[HE]): 

V^pTr^) = p Tr{^ . $ (2) • > 0, 

where the positive semi-definiteness comes from the fact that p > and 7y and $ are 
positive definite. Denoting by dT 1 the vector surface element on dT pointing from T 
outwards, the Stokes theorem applied to the above identity leads to 

/ p Tr * ti dV = f p 7^ Tr{^ . $ . j-i . $ }^ ( 27 ) 

as long as pTr^jdY* = at T^, and the integrand on the left is continuous up to dT. 
Now, defining / = + / (2) , one has [TJ1 1201 E] 

1 



Tr^ 



X(i)X {2) 



Xf + 2X\ Ar + I A| 4 + Y + Im(AA) 



(28) 



Taking the limits of the values of X, Y and A on the axis and at infinity -with decays 
given by (P)- (fTU|) . but independently of the values of M, J, and Q- of the configuration 
as were computed by Carter 20 (see also |T21 El), one first obtains that Tr\& and 
Tr tyj i dT l \r ± are regular (C 1 ) on the axis 20J, so that pTr \J/ jdT' J |x ± = 0, and secondly, 
that p Tr ^dY^x^ also vanishes [20]. As a result, pTj) holds and the only contribution 
to the integral on the left hand side can come from Y^. On the other hand, since the 
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numerator of Tr \1/ is quadratic in the differences of the variables, X, etc, and $|t^ =0 
by assumption, one infers dTr\l/|-r s = 0, and hence p Tr ^dT^gr = 0. 

Therefore, due to the positive semi-definiteness of the integrand on the right hand 
side, (J2ZJ) implies that J_ = in T (and thus in W E ). Since d^ = $ ( i) • J • $(2) -1 by 
construction, \1> is constant all over T, and thus \1/ = because \I/|t e = by assumption 
(J2HJ). This ends the proof showing that given Dirichlet data on Ts, and correspondingly 
on S, i.e. {^|s,A^| s }, the solution A^) of the Ernst equations in the exterior 
region (W E ,g E ) is unique. The correspondence between conjugate solutions (£^>, A^) 
and (£g, A|*) (see jHj) leads to the same result for k = £. The main result found in this 
section can be stated as follows: 

Proposition 1 Let (W E , g E , !F, X) define a simply connected stationary and axially 
symmetric electrovacuum spacetime with a connected boundary preserving the 
symmetries. If W E is asymptotically flat and the problem is given Dirichlet boundary 
data on E for the Ernst and electromagnetic potentials, i.e. A|s}, then the solution 

(£, A) is unique. ■ 



5. Fixing the Dirichlet data 



The purpose now is to show how the full set of Cauchy data, i.e. taking into account 
{d£| s , dA| s }(co>, A) on top of the family of Dirichlet data A| s }(cu, A), fixes the 
values of u, A, provided that the solution exists. 

The proof presented here makes use of the divergence free fields 3 (OHj) choosing 
k = £. Decomposing J" ^ on a basis of su(2, 1) one can obtain eight conserved real 
1-forms. For convenience, the basis used here corresponds to that presented in [T^j, so 
that the eight forms will be presented as a set of four real 1-forms and two complex 
1-forms. The representation of an element Z G su(2, 1) = ee + 5d + As + np + ah + (3a, 
where four elements of the basis of the algebra are denoted as e, d, s, p and the remaining 
four are encoded in h and a by "complexifying" them, is chosen as ^H] 

( -V2a e \ 



6 {Z) 



2" 3 
TV 



2 \ 



(29) 



i^/2 a 

where a and (3 are complex contants and the rest are real. Nevertheless, one cannot 
yet revert to the decomposition of the divergence-free currents in ^Hj, because the 
construction of the currents satisfying (|2"2|) is performed in |UU in a different way. In 
fact, whereas the vacuum case lies on su(l, 1) when using (|2^ j) - (j2^j) . the construction 
presented in [TH] leads to an isomorphic sl(2, 1R) in the vacuum case. Denoting by J~ BM 
the matrix representation of the divergence-free currents as constructed in JH] for k = £, 
and defining 

( 1 -i\ 



a 



\ 










i 



J 
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the isomorphism between the two constructions can be found to be given explicitly by 

J^ = \o" ■ J v o. (30) 

Using (j2HJ) »7bm is decomposed (see also jT!|| pp.798) in terms of a basis of four real 
1-forms J7" e , J7" d , J7" s , 3 V plus two complex ones, J~ h and J7 a , which corresponds to the 
elements of the basis of the algebra su(2, 1) listed above, as 

Jbm = \j e e + J d d + J s s + ~J>p + ~ij h h + \j a a. 

The explicit expressions, which can be obtained from the above construction (|22j) -()23 j) 
by making use of (j3*Uj) . are given in Appendix B| Using these 1-forms, one defines on 



(r, 7) a (complex) 1-form depending on a (£, A) configuration plus eight real parameters 
as 

W^(£, A, c e , Cd, c s , Cp, c/jj , c/j 2 , c ai , c a2 ) = 

CeJ" + C d J d + C S J S + C p J P + C hl J h + C h2 J h + C ai J a + C a2 J a , 

where clearly all J = J{£, A), which by construction and (J2EJ) satisfies 

w V*(pWi) = 0. (31) 

The only non- vanishing surface integrals at spatial infinity of the above 1-forms 
are given by (see e.g. [TH] ) 

/ p J d i dT* = AM, f p J h i dT= [ p J\ dT* = 2Q, 

where M(£,A) and Q(£,A) relative to a given solution (£,A) correspond to the mass 
and the electric charge of the configuration, respectively, so that at spacelike infinity 
one has (J9~|) - ()l()j) . As a consequence of (|31j) we have 

- / p W t dT = [ p WAT = 4Mc d + 2Q (c hl + c ai + c h2 + c a2 ) . (32) 

Now, let us assume two exterior solutions for the same interior exist (5m, Am), 
(5(2), A(2)), with corresponding Mm, Q(i) and M( 2 ), Q(2) respectively, such that their 
Cauchy boundary data differs by A and to as given by the relations (JT2J), fl~3|) . (}20|) . and 

A(2)|t e = A(i)|x E + A, 5(2)|x E = £(i)|t e — AA + iuj — 2AA| Te , 
dA( 2 )|T E = dA(!) | Ts , d£( 2 )|T E = d5(i)|x E — 2AdA| Ts . 

From these relations, the following equality holds on T^, 

W(£m,Am;c e , ...,c s )| Ts = W(5(2),A(2);c e , ...,c s )|x E) (33) 
for a set of eight certain relations 

c e = c e (c e , . . . , c s ; cj, A), . . . , c s = c s (c e , . . . , c s ; o>, A), (34) 
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which are given explicitly in Appendix B The integration of (|33|) over using (|32j) 
leads to 

4M(i)Q + 2<5(i) (c hl + c 0l + c h2 + c a2 ) 

= 4M( 2 )Q + 2<5 (2 ) (c/ji + C„i + c h2 +Ca 2 ) , (35) 

which has to hold for arbitrary choices of the eight parameters {c e ,...,c s }. A 
straightforward calculation, see |Appendix B , leads to the fact that 



M ( i) + M (2 ) ^0^A = c; = 0, 

(so that in fact M(i) = M( 2 ), Q(i) = Q(2) necessarily). It must be stressed here that 
the case M(i) + M( 2 ) = does not lead to the same result in general, see Appendix B 



Nevertheless, and since the mass for physically motivated solutions is positive, this result 
basically implies uniqueness of the Cauchy data for which the exterior field exists: the 
freedom in the Cauchy and Dirichlet data generated by A and to is fixed. To be more 
precise, we have proven the following: 

Proposition 2 Let us consider (W E , g E , J- , X) as in Proposition^ with a family of 
Dirichlet boundary data depending on three real parameters, i.e. A|e}(o;, A) ; so 

that it is determined up to transformations of the form (Q3j) and A2U\) . Consider now the 
completion to Cauchy boundary data by adding the data {d£|s, dA|s}(w, A) determined 
up to transformations of the form dA'| s = dA| s and $H\) . If within the (tu, A) -family 
of Cauchy data sets there is a Cauchy data set for which a solution with M > exists, 
then there is no solution with M > for any other Cauchy data set within the family. 



6. Conclusions 

The combination of propositions ^ and |21 with the boundary data that a given interior 
infers onto the exterior problem, as shown in Sectional leads us to: 

Theorem 1 Let (W 1 ,^ 1 ) be a given stationary and axially symmetric region of 
spacetime containing an electromagnetic field, and bounded by a symmetry-preserving 
and simply connected hypersurface S ; spatially homeomorphic to a two-sphere, 
and containing no surface charge currents. If (W 1 ,^ 1 ) can be matched across 
£ to an asymptotically flat (with M > 0) stationary and axially symmetric 
electrovacuum Einstein-Maxwell field (W E ,g E ,^F), preserving the symmetry, and where 
the identification on E has been prescribed, then (yV E ,g E ,^F) is unique within the set of 
positive mass solutions. ■ 

As a first remark, the same result holds if £ is allowed to contain surface electromagnetic 
charge currents, as long as these currents are taken as given data. Secondly, although 
for physically motivated situations the restriction of the uniqueness result to the set of 
positive mass solutions is irrelevant, it may be desirable to have a complete uniqueness 
result. This might be expected to follow by using existence considerations, and it is 
currently under study. 
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All in all, this final result can be stated in other words as follows: under 
natural assumptions on X -preservation of the symmetries and the prescription of the 
identification of both sides of £ [2] , the exterior electrovacuum field generated by an 
isolated distribution of charged matter in axial rotation and in equilibrium is unique, 
provided that it exists. 

Some comments on the definition of asymptotic flatness are in order here. The 
asymptotic flatness has explicitly entered the problem only through the use of the 
asymptotic conditions (PJ)- ([TUj) . In principle, one could have included NUT and 
monopole charges in the expressions for U and A at infinity. In fact, it can be shown 
that the inclusion of NUT and monopole charges does not change the present results, 
provided we can use Weyl coordinates globally on (W E , g E ) (see (J2J)). Nevertheless, the 
existence of global Weyl coordinates has been only ensured when (W E ,g E ) admits an 
asymptotically flat four-end [Hj (see e.g. [21] for definitions). Therefore, this is the 
definition of asymptotic flatness that has been implicitly used all throughout the paper, 
which implies the vanishing of the NUT and monopole charges as mentioned in 
Section |U 

Looking at the procedure used here, several generalisations to include different 
fields on the exterior region could be attempted quite straightforwardly. The starting 
point would be the possibility of using the Weyl coordinates (orthogonal transitivity 
and a trace-free energy- momentum tensor [5 ). Then, one should consider the boundary 
conditions for the fields on the matching hypersurface, to be combined to those inferred 
by the spacetime matching conditions, according to Proposition EJ in order to find 
the form of the Cauchy boundary data for the exterior problem. Of course, the main 
ingredient would consist of picking up fields that can be described by a-models [10], 
and use then the Mazur construction for the Dirichlet problem, and the divergence-free 
currents for the fixing of the Dirichlet data with the whole set of Cauchy data. 
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Appendix A. The explicit matching 

This section is devoted to explicitly present the general matching of two stationary and 
axially symmetric spacetimes (with boundary), (W 1 ,^ 1 , S 1 ) and (W E , g E ,£ E ), across 
a matching hypersurface S (= X 1 = S E ) preserving the stationarity and the axial 
symmetry. The latter implies that £ is timelike. Furthermore, we impose on (W E ,g E ) 
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the following two requirements: (i) the G 2 on T 2 group of isometries acts orthogonally 
transitively (OT), and (ii) the trace of the Einstein tensor is zero. These two conditions 
simply account for having the possibility of simplifying the problem by using the Weyl 
coordinates (at least locally) (see j^j). 

Apart from preserving the symmetries, the only implicit assumption we are making 
on E is that it intersects points at the axis, so that the axis intersects both W 1 and W E , 
which is ensured if E is spatially homeomorphic to a two-sphere. Conditions (i) and (ii) 
are satisfied when (W E ,g E ) contains a stationary and axisymmetric electrovacuum field 

u 

The equations for the matching with a vacuum (W E , g E ) spacetime where presented 
in an explicit manner in [7j. The reader is referred to that work (and also to 0122] for 
the OT (W 1 , g l ) case) for further details. 

There exists a coordinate system {T, $, r, (} in which the line-element for g l reads 

[23 El 

ds 2 = -e 2V (dT + 5d$ + W 2 dC) 2 + e~ 2V [e 2h (dr 2 + d( 2 ) + a 2 (d$ + W 3 d() 2 ] , (A.l) 

where V, B, h, W 2 , W 3 and a are functions of r and (, and ff 1 = d$ is the (intrinsically 
defined) axial Killing vector. The vector £ : = &t is a stationary Killing vector. By 
assumption, (W E ,g E ) admits Weyl coordinates {t, <f>, p, z}, in which its line-element is 
given by (j2J), the intrinsically defined axial Killing vector is given by ff E = d^, and we 
take £ E = dt. As mentioned in Section 121 it is straightforward to show that the most 
general matching hypersurface preserving stationarity and axial symmetry E, as seen 
form the (W E ,g E ) side, is given in parametric form as [2} II 1| IT^j 

£ E : {t = r, <p = <p, p = p(/i), z = z(p)}, 

where {r,ip,p} parametrise E. At this point, the matching procedure introduces two 
degrees of freedom, a product of the identification of dt with any linear combination 
a(dx + bd§) on E [21 HU Ej. Of course, this is of relevance if d t is intrinsically defined, 
as it happens, for instance, for asymptotically flat (W E ,g E ) j2J. Despite that, one can 
proceed with the matching avoiding the explicit appearance of a, b by absorbing them 
in (W 1 , (7 1 ), changing coordinates in (W l ,g l ) to a new set of "primed" coordinates, with 
corresponding functions V, B', h', W 2 , W% and a' in (jA.lj) . so that 8t> = a(dr + bd<$>) 
while keeping 9$/ = <9$ (see [7j for the explicit relations between primed and the 
original non-primed quantities). Nevertheless, if the interior geometry (jA.lj) was given 
beforehand, these relations involving the parameters a, b have to be taken into account. 
Having this in mind, and dropping the primes, the most general parametric form of E 
as seen from (W 1 , g l ) reads then [2] 

E 1 : {T = r + f T (p), $ = <p + f (p), r = r(p), ( = ((//)}, 

where f T (jj) = ((BW 3 - W 2 )(p) and U(p) = -(W 3 (p). 

Two basis for the tangent spaces to E, from the E 1 and E E sides, {e l a } and {e E } 
a : 1,2,3 respectively, which are to be identified eventually as {e a } in the final matched 
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spacetime (21], read explicitly 

e{ = £%, el = ff%, e\ = rd r + ( [(BW 3 - W 2 )d T - W 3 d* + d c ] | s , 
ef = f E |s, e 2 E = rf | E , e 3 E = pd p + i<9 2 | 2 . 

In principle, any choice of normal vector n is suitable for computing the matching 
conditions, but it is fixed, by convention, so that it points (W 1 , g 1 , S 1 ) inwards and 
(W E ,g E ,S E ) outwards, and by convenience, that its norm is the same as that of e 3 . 
On the S E side it reads n E = e 2 ( k ~ u \—z dp + pdz), while on the S 1 side, this is 
n 1 = e e 2<yh ~ v \— C, dr + fd(), where e is a sign. The equality of the first fundamental 
forms is equivalent to the following four equations, the so-called preliminary junction 
conditions, 

U\v = V\v, A\ E = B\v, (A.2) 

p|s = a|s, (A. 3) 

e 2fc | s (p 2 + i 2 )=e 2fe | E (r 2 + C 2 ). (A.4) 

Since £ is timelike the remaining matching conditions simply consist of equating the 
second fundamental forms, K ab = e^e^V ' a np, which are symmetric [21], as computed 
from E E and S 1 . Labelling by A, B = 1, 2 the components corresponding to the orbits, 
the set of six equations K^ b = K l ab , after using (|A.2J) - (|A.4J) . can be expressed and grouped 
as follows: 

K1 B = K\ B : n*(U)\x = n\V)\ s , n E (A)| E = n\B)\ s , (A.5) 
z = -f?(a)| E (A.6) 
K% = K\ Z : d r (BW 3 -W 2 )\ s = 0, d r W 3 \^ = (A.7) 

K! 3 = K 33 : e 2k (p 2 + z 2 )n E (k) + e 2 ^ (pz - pz) | s 

= e 2h (r 2 + C, 2 )n\h) + e 2 ^ (f( - r() | E , (A.8) 

where flA.5|) has been used in (jA.8|) . 

Using all the equations one finds first that p 2 + z 2 = (a 2 r + a^^ij 2 + C 2 ), meaning 
that a 2 r + a 2 Q cannot vanish anywhere on S, and secondly, that (jA.4|) and (jA.8|) can be 
substituted by 

1 



fc|s = h- - log(a 2 r + a 2 f ; 



s 



n*{k)U=n\h)-e ^l ■ "f' 

ay + 



E 



(A.9) 
(A.10) 



respectively. A straightforward calculation shows that, once the rest of the equations 
hold, equations (jA.10|) and the derivative of (jA.9|) along p are equivalent to 

n Ba n B ^| E = n^n^S^, n* a efS* p \x = n^efS^ (A.ll) 

where S a p denotes the Einstein tensor. On the other hand, the equations in (IA.7|) are 
equivalent to [7j 

*(? Ari' A dOls = = *(e A rf A drj l )\^ 
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(see also |TT|). and, in turn, equivalent to [7j (up to additive contants, but these vanish 
because £ intersects the axis) 

n la e I 1 f, S l af) \ s = 0, n la e I /S l a/3 \x = 0. (A.12) 

Since n la el P Sl p = = n la e 2 S^p are identically satisfied due to the OT structure of 
the (W E , g E ) region, the two equations in (jA.12|) together with those in (|A.11|) constitute 
the set of the four Israel conditions. 

To sum up, the whole set of matching conditions is formed by the ten equations 
in (lA~2l) . (TPL (lA~5|) . (TPL (jOJ), (lA~Toi and (|A~T2jl . taking into account that up to 
an additive constant in k\s, equations ()A.9j) - (jA.10j) can be substituted by (jA.llj) . If 
required, that constant should then be fixed by considering (|A.9|) . 

In principle, if no conditions are imposed on the matter content in (W E ,g E ), 
equations (jA.9|) - (jA.10j) . or equivalently (jA.llj) . provide conditions on £ for the function 
k, given the interior geometry. Conversely, if one imposes conditions on the matter 
content in (W E ,g E ), i.e. on S E g, so that the function k has to satisfy equations 
involving the rest of the functions, equations (|A.9|) - (jA.10|) may represent conditions 
on the (W 1 , g 1 ) side only. In particular, this becomes clear when (W E ,g E ) is vacuum, 
because (jA.ll|) only involve quantities in (W 1 ,*? 1 ), and therefore ()A.9|) - ()A.10|) consist, 
in fact, of conditions on S 1 (apart from providing the value of the additional constant 
in k). 

Although not as explicitly apparent as in the vacuum case, it is straightforward 
to show that the same holds true when (W E ,g E ) is an electrovacuum solution whose 
electromagnetic field extends continuously across £ onto (W\ g 1 )- As mentioned in 
Section 121 k is determined, up to an additive constant, by quadratures in terms of the 
rest of the metric functions, which satisfy the Ernst equations. Defining T = — dN + iw, 
where this and the following quantities are referred to k = £, the quadratures for k in 
((21) can be cast as (see 0E!) 

~ p K* = (r„f z + T Z T P ) + i (A, P A, Z + A,A, P 

The explicit expressions in terms of U and A are then obtained by using and a 

straightforward calculation on S leads to 



kk= ^—(Ap-Bz) 



ti*(k)\x= ^—(Az + Bp) 
p A + z £ 



where we have defined 

A = U 2 + n E (U) 2 + [A 2 + n E (A) 2 ) - e~ 2U (AA - n E (A)n E (A 



e 4U 



4/9 

B = 2Un E (U) - ——^An E (A) - e~ 2U (n E (A)A + An E (A) 



2p 2 

The continuity of the electromagnetic field across S implies relations for A|s and 
n E (A)| s in terms of (I) quantities. Using these relations together with (|A.2|) . (jA.3l) . 
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(|A.5J) and (|A.6J) . and n E (/c)|s can be obtained in terms of quantities of (W 1 , g 1 ), and 
hence obtain expressions for (jA.10|) and the derivative of (jA.9|) along fi as conditions 
for quantities on the (W 1 , y 1 ) side only. Equivalently, equations (jA.llj) . within the set 
of matching conditions, provide conditions on (W 1 , p 1 ) only. 

On the other hand, of course, it must be stressed that conditions on the matter 
content in (W E ,g E ) only translate to conditions on the (W 1 ,*? 1 ) geometry through the 
four Israel conditions (jA.ll)) - (jA.12j) . and thus, given an interior geometry satisfying 
these conditions, the rest of conditions involve the form of the matching hypersurface 
at the b side, E E , plus boundary data (of Cauchy type) for U and A on E E given the 
interior geometry. 

All the above can be summarised as follows: 

Proposition3 Let (W E ,g E ) and (W\ g 1 ) be two stationary and axially symmetric 
spacetimes which are to be matched across their respective symmetry-preserving 
boundaries E E and E 1 , whose identification is denoted as E(= E E = E 1 ). Let us also 
assume that (W E ,g E ) contains a stationary and axially symmetric electromagnetic field 
and satisfies the Einstein-Maxwell equations without sources. The whole set of matching 
conditions can be reorganised as follows: 

(i) Conditions on E 1 , given by hA.ll\) . \A.12i) . which form, in principle, an 
overdetermined system for {r(/z), C(/-0}- If « solution exists the matching is possible 
(locally) and the hypersurface E 1 is determined generically. 

(ii) Equations defining E E ; given by \A.$) . KA.fy) . 

(Hi) Boundary conditions for the electrovacuum problem in (W E ,g E ) ; given by iA.fl) . 
(L4..5)) . If required, the additive constant in k is determined by using ( L4..9)) . ■ 

Remark 3.1 The explicit expressions apply for the definitions and metric functions as 
introduced above in this Appendix, but can be easily reexpressed in an intrinsic manner. 



Appendix B. The divergence-free forms 



The explicit expressions of the divergence free (in (W, g)) 1-forms used for the 
decomposition of the su(2, l)-valued 1-form (for either k) are given as follows, both 
in terms of 8 and A and in terms of the sometimes more convenient N, A, fl and w (for 

dn (HD): 



2N 2 
1 



N 2 



w. 



2N 2 
1 

3i 



AN 2 
3 



2(AdA - AdA) - d(£ - S) 
d(£S) + AAd(£ + £)- (AdA - AdA)(^ - 8) 
VLw + N(dN + d(AA)) 

AAd(^ -£)- (AdA - AdA) (8 + 8) 



2N 2 



AAw - iN(AdA - AdA) 
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J p 



ja 



2N 2 
1 

1 

2N 2 
1 

2N 2 



2££(AdA - AdA) + (£ + 2AA)£d£ - (£ + 2AA)£d£ 



iV 2 - ft 2 - A 2 A 2 )™ - 2Nft dN - 4Re 



i(N + AA)) NAdA] } 



Ad(£ -£) + {£ + £)dA + 2A 2 dA 



N 2 



iAiu - ATdA 



£{£ + £)dA + 2A 2 (£dA - Ad£) - Ad(££) 



— {A^ 2 dA - A^AdA" - 2NA 2 dA - (fi + zAA)(Au> + iA^dA)} . 



AT2 

The vacuum case is recovered by A = (so that dQ = w), and thus only three 1- 
forms survive, namely J e , J d and J k . Their expressions in the = £, and 

after the substitution N = -e 2U , read J e = e~ 4U dn, J d = 2dU + e~ 4U ndn, and 
ST P = — 4QdU + (1 — e~ 4U Q 2 )dQ, which were those used in constructing the divergence- 
free 1-form in [2|. 

With the above expressions one is ready to compute the relations (j^jj) by using the 
differences in the Cauchy data of the two solutions in terms of A, N, Q and w l|12 J) - p9|) : 



A (2)|t e 


= A (i)lx E + A, 


dA(2)|r E 


= dA(i)| TE , 


A^( 2 )|t e 


= ^(i)It e , 


dA^(2) t e 


= d%)| Ts , 


^(2)|t e 


= ^(i)|x E +w + i(AA(i) - AA (1) )| Te , 


W (2)W S 


= W(i)| Ts - 



Indeed, is recovered with 

c p = Cp, c a i = c i + i2c p A, c a 2 = c a 2 — i2c p \ 

Cd = c d + c al X + c a2 A + 2c p uj, c s = c s + 3i(c a i\ - c a2 A) - 6 c p AA, 

41 = c hl + (c d + ic s )A - c ai (AA + iui) + 2A 2 (c a2 - ic p X) + 2c p A^, 

4 2 = c/i2 + ( c d - ^c s )A - c a2 (AA - Lui) + 2A 2 (c ai + ic p X) + 2c p Xu, 

c e = c e + i(c h2 X - c hl X) + XX(ic ai X - ic a2 X - c p XX + c s ) - (c d + c ai A + c a2 A + c p u)u. 

These expressions are then introduced into ()35)). which has to hold for arbitrary 
{c e , . . . , c s }. Putting Ch x = 1 and the rest zero (say, the c^-equation: in what follows, 
the analogous procedure will be referred to as c x -equations) one immediately obtains 

Q(i) = Q{2), 

whereas the equation involving c s implies 
Q( 2 )(A- A) = 0. 

Using these two relations on a suitable combination of the equations involving c ai and 
Cd, explicitly A[A(crf-equation)— 2 (c ai -equation)], one obtains 

2AA(M (1) + M (2) ) = Q {2) ui(X + A), 

whose real part reads simply 

2AA(M ( i) +M (2) ) = 0. 
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The same procedure applied to another combination of the equations involving q and 
c p , explicitly cj(Q-equation)— 2 (c p -equation), leads to 

2w(M ( i) +M (2) ) = -wQ(a)(A + A). 

From the last two equations, if Mm + Mm 7^ 0, then A = u = 0, and thus the main 
result in Section El follows. As expected, once A = u = it is immediate to see that 
then M (1) = M (2) . 

It must be stressed here that the case Mm + M( 2 ) = does not lead to the same 
result: by setting M( 2 ) = — Mm 7^ 0, it is immediate to find that uj = 0, A — A = (so 
that only the real part of A survives) and finally, A = 2Mm/Qm exhausting equation 
(HU). 
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